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ABSTRACT 



Very high energy (VHE) 7-ray emission from the flat spectrum radio quasar 
(FSRQ) PKS 1222+21 (4C 21.35, ^=0.432) was detected with the MAGIC 
Cherenkov telescopes during a short observation (~0.5hr) performed on 2010 
June 17. The MAGIC detection coincides with high energy MeV/GeV 7-ray ac- 
tivity measured by the Large Area Telescope (LAT) on board the Fermi satellite. 
The VHE spectrum measured by MAGIC extends from about 70 GeV up to at 
least 400 GeV and can be well described by a power law dN/dE oc E"'" with 
a photon index F = 3.75 ± 0.27^^^^ ± 0.2,y,t. The averaged integral flux above 
100 GeV is (4.6 ± 0.5) x lO'^Ocm-^s"^ (~ 1 Crab Nebula flux). The VHE flux 
measured by MAGIC varies significantly within the 30 minutes exposure imply- 
ing a flux doubling time of about 10 minutes. The VHE and MeV/GeV spectra, 
corrected for the absorption by the extragalactic background light (EBL), can be 
described by a single power law with photon index 2.72 ± 0.34 between 3 GeV 
and 400 GeV, and is consistent with emission belonging to a single component 
in the jet. The absence of a spectral cutoff constrains the 7-ray emission region 
to lie outside the broad-line region, which would otherwise absorb the VHE 7- 
rays. Together with the detected fast variabihty, this challenges present emission 
models from jets in FSRQs. Moreover, the combined Fermi /LAT and MAGIC 
spectral data yield constraints on the density of the EBL in the UV-optical to 
near-infrared range that are compatible with recent models. 

Subject headings: cosmic background radiation — galaxies: active — galaxies: jets — 
gamma rays: galaxies — quasars: individual (PKS 1222-1-21) 
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Introduction 



High- luminosity active galactic nuclei (AGNs) hosting powerful relativistic jets are 
characterized by strong nonthermal emission extending across the entire electromagnetic 
spectrum, from radio up to 7-rays. More than 40 AGNs have been detected in the very high 
energy (VHE) domain {E > 100 GeV) by ground based Cherenkov telescopei^]. The great 
majority of them are BL Lac objects, whil e only two are classified as fiat spectrum radio 



quasars (FS^ 


IQs): ] 


Albert et al. 


2008a 



'z = .36, IWagner fc Behera 



Aleksic et al 



2OI0I) and 3C279 {z = 0.536, 



2OIII ). the most distant VHE source detected up to now. 



FSRQs display luminous, broad emission lines often accompanied by a "big blue bump" 
in the optical-UV region, associated with the direct emission from the accretion disk. VHE 
emission from FSRQs may therefore be affected by internal absorpt ion from the dense 



UV-o ptical radiation reprocessed in the Broad Line Region (BLR) ( iDonea fc Protheroe 
20031 ). Distant VHE quasars offer the possibility to probe the Extragalactic Background 
Light (EBL ), the integrated stella r and dust emission through cosmic history, in the range 



0.1-10 fim (IHauser fc Dwek 



2OO1I). 



The MAGIC detection flMariotti et al. 
z = 0.432 



Osterbrock fc Pogge 



2010h of the FSRQ PKS 1222+21 (4C 21.35, 



19871 ) makes it the second most distant obie ct with know n 

Jl ' 1 

redshift (after 3C 279) detected at VHE^I. PKS 1222+21 is a 7-ray blazar fAbdo et al. 



2010bl ) with a relatively hard spectrum in the GeV range and has been included in the list 



Neronov et al. 



(|2010[ ). It is characterized by highly 



of >100 GeV emitters in the analysis of 

^For an updated list refer to 
http : / / www. mppmu . mpg.de / ~rwagner / sources / 

^ The redshift measu rement [z = 0.444) of the VHE BL Lac 3C 66A has large uncertain- 



http: / /tevcat. uchicago.edu/ 



or 



ties ( Bramel et al 



2OO5I). 
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superluminal jet knots with apparent velocity up to 21 c flLister et al.ll2009[ ) 



Uppe r limits on the VH E emission of PKS 1222+21 have been previously derived by 



Whipple (IKerrick et al. 



19951 ) at the level of 12 x 10-^2^^-23-1 at E> 300 GeV. We report 
here on the MAGIC discovery of this source, during a phase of high activity in 7-rays 
announced by the Fermi /LAT collaboration. We discuss its implications for the EBL 



studies and the blazar physics. 



Observations 



MAGIC consists of two 17 m diameter Imaging Atmospheric Cherenkov Telescopes 
(lACT) located at the Roque de los Muchachos, Canary Island of La Palma (28°46' N 



17°53'W), at the height of 2200m a.s.l. The stereo observation s provide a sensitivity 



0.8% of the Crab nebula flux at E> 250 GeV fjCohn et al. 



20091 ). 



^lof 



PKS 1222+21 was observed by MAGIC from May 3 to June 19 2010 (MJD 55319 
to MJD 55366) for a total of ~14.3 h. The observations started as a part of a Ta rget of 



Oppo rtunity program triggered by an increase of the flux in the Fermi passband (.Donato 



20101). In this letter we report the results obtained from the observation of June 17 (MJD 



ose coincidence wit h the brightest 



55364), when the source was detected by MAGIC in c. 

flare observed by Fermi Large Area Telescope (LAT) ( iTanaka et al.ll201ll ). Results from 
the multi- wavelength campaign covering all 2010 observations will be published elsewhere. 
Nevertheless a preliminary analysis does not provide any high-significant detection with 
MAGIC in any other day during the campaign. 



^Sensitivity is defined here as the minimal integral flux to reach 5a signal in 50 h of observations. 
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On June 17, 21:50 UT, PKS 1222+21 was observed with the MAGIC telescopes for 
~0.5h (MJD 55364.908 to MJD 55364.931), in the so-called wobble mode. The data were 
taken at zenith angles between 26° and 35°. The light conditions during the observations 
correspond to moderate moon light leading to a higher noise level in the data. A cleaning 
level higher than the standard one was therefore applied to remove signals from night sky 
background noise. Stereoscopic eve nts, triggered by both MAGIC telescopes, were analyzed 



in the M ARS analysis frame work (IMoralejo et al 



found in 



Aleksic et al. 



20091 ). Details on the analysis can be 



(|2010[ ) whereas the performance of the MAGIC telescope stereo 



system will be discussed in detail in a forthcoming paper. 



Results 



The strength of the signal was evaluated applying standard cuts to the PKS 1222+21 
data sample, corresponding to an energy threshold of ~ 70 GeV as determined by Monte 
Carlo events, assuming a soft spectrum with a photon index of F = 3.5. The fP' distribution 
(squared angular distance between the true and reconstructed source position) of the signal 
coming from the region of PKS 1222+21 yields an excess of 190 - r-like events (67 /min.), 



corresponding to a statistical significance of 10.2 a using eq. 17 in 



Li fc Mai fll983h . 



3.1. VHE spectrum 



The differential energy spec trum of PKS 1222+ 21 was reconstructed using the 



'Tikhonov" unfolding algorithm (lAlbert et al. 



20071 ). to take into account the finite energy 



resolution of the instrument and the biases in the energy reconstruction. The energy 
spectrum, shown in Fig. [H extends up to at least 400 GeV and is well-described by a simple 
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Fig. 1. — Differential energy spectrum of PKS 1222+21 as measured by MAGIC on 2010 
June 17. Differential fluxes are shown as black points, upper limits (95% CL.) as black 
arrows. The black line is the best fit to a power law. The grey shaded area represents 
the systematic uncertainties of the analysis. The absorp tion corrected spectrum and upper 



limits using the EBL model by iDominguez et al. 



( I2OIII ) are shown by the blue squares and 



arrows; the dashed blue line is the best fit power law. The blue-striped ar ea illustrates the 



uncert a inties due to d i fferen c es in the EBL models cit ed i n the text, by 



mm; 



Gilmore et al 



( l2009h : 



Franceschini et al 



f l2008h and 



Albert et al 



Kneiske fc Dole 



f!2008ah 
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power law of the form: 



— -N ( ^ 
dE ~ V 200 GeV 



with a photon index F = 3.75 ± 0.27^4^1 =t 0.2gygt and a normahzation constant 
at 200 GeV of A^aoo = (7.8 ± 1.2,,,^ ± 3.5,^,0 x IQ-^^cmrh-^TeV-^ , yielding an 
integral flux (4.6 ± 0.5) x 10"^° cm-^g-i 1 Crab Nebula flux) at E> 100 GeV and 
(9.0 ± 3.6) X 10-^2 cni-2s-i (7% of the Crab Nebula flux) at E> 300 GeV, at the same 
level of Whipple upper limit (sec.l). For energies higher than 400 GeV no significant excess 
was measured. The upper limits corresponding to 95% confidence level (C.L.) are shown 
in Fig. [TJ The systematic uncertainty of the analysis (studied by using different cuts and 
different unfolding algorithms) is shown by the grey area. 

We studied the effect of the VHE 7-ray absorption due to pair-production with 
low energy ph otons of the EBL by usi n g different state-o f -the-a r t EBL models, namely 



the models by 



Dominguez et al. 



Franceschini et al. 



(1201lh : 



Kneiske fc Dole 



(I2OI0I) 



Gilmore et al. 



(120091): 



(120081 ) and the "max high UV" EBL model described in lAlbert et al. 



(j2008al ). For each of the EBL models the optical depth corresponding to the measured VHE 
7-ray energy intervals was computed and the differential fluxes were corrected accordingly 
to obtain the de-absorbed ( or int rinsic) spectrum. The spectrum deabsorbed with the EBL 



model of 



Dominguez et al. 



(I2OIII ) , shown by the blue squares in Fig. [U are well fitted by 
a power law with an intrinsic photon index of Fi^tr = 2.72 ± 0.34 between 70 GeV and 
400 GeV. Uncertainties caused by the differences between the EBL models are represented 
in Fig. [T] by the blue-striped area. The corresponding spread is smaller than the systematic 
uncertainties of the MAGIC data analysis. 

We investigated the possible presence of a high energy cut off in the VHE range by 
fitting power laws with different photon indexes and d ifferent values for the cut off. The 
method adopted is the difference method (see, e.g. 



Lampton et al. 



19761 ). With the 
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available statistics, at the 95% C.L. we cannot exclude the presence of a cut off above 130 
GeV for a photon index 2.4 (the lowest possible value compatible with fit uncertainties 
and with Fermi /LAT data, see sec. 3.3) or above 180 GeV for a photon index 2.7. The 
confidence interval is not bounded on the high energy side, i.e. a fit without a cutoff is fully 
compatible with the data. Further observations with higher statistics are needed to better 
constrain the location of a possible steepening in the form of a cut off or spectral break. 



3.2. Light curve 

Despite the short observation time, the strength of the signal allows us to perform a 
variability study of the measured integral fluxes above 100 GeV. The light curve binned in 
6 minutes long intervals is shown in Fig. [2] and reveals clear flux variations. The constancy 
hypothesis {x^/NDF = 28.3/4) is rejected with high confidence (probability < 1.1 x 10~^). 
The fluxes of background events surviving the 7/hadron selection cuts are compatible with 
being constant and hence we can exclude a variation of the instrument performance during 
the observation. 

To quantify the variability time scale we performed an exponential flt (solid black 
line in Fig. [2]). A linear flt is also acceptable but does not allow us to deflne a time 
scale unambiguously. For the exponential flt the doubling time of the flare is estimated 
as 8.6]!^o 9 '^i'^^t^s. The derived timescale corresponds to the fastest time variation ever 



observed in a FSRQ in the VHE range and in any other energy range (iFoschini et al. 
and i s among the shortest timescales measured on TeV emitting sources ( jAbramowski et al. 
2010k 
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Time (UT) 



Fig. 2.— PKS 1222+21 light curve above 100 GeV, in 6 minutes bins (black filled circles). 
The observation was carried out on MJD 55364. The black solid hne is a fit with an ex- 
ponential function and the black dotted line a fit with a linear function. The grey open 
squares denote the fluxes from the background events and the grey dashed line is a fit with 
a constant function to these points. 
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3.3. The HE - VHE SED 



In the high energy (HE) MeV/GeV energy range measured by Fermi /LAT the source 



showed a significant 



flare lasting ~3 days, with a fiux peak on 2010 June 18 (MJD 



55365) (ITanaka et al. 



1201 ll ). A dedicated analysis found that the 1/2 h MAGIC observation 
fell within a gap in the LAT exposure, thus we analyzed a period of 2.5 h (MJD 55364.867 
to 55364.973), encor npassing the MAGIC observation. The LAT analysis for this time bin 



was performed as in iTanaka et al. 



(120111 ). where details can be found. It results in an 
integral flux (6.5 ± 1.9) x 10~® cm~^ s~^ at energy E> 100 MeV. The observation in such a 
short time does not provide any detection with Fermi /LAT at E>2GeV. Two Fermi /LAT 
spectral points up to 2 GeV together with an upper limit at the 95% C.L. in the range 
2 — 6.3 GeV are combined with the MAGIC data in the Spectral Energy Distribution (SED) 
shown in Fig. [3l 

The flgure also shows bow ties representing uncertainties associated with the spectral 
flts. The Fermi /LAT spectrum is best described by a single power-law with index of 
1.95 ± 0.21. In the case of MAGIC data the bow tie refers to the "intrinsic" source 
spectrum, i.e. to the observed spectrum corrected for EBL absorption, described in sec. 3.1. 
An extrapolation of the intrinsic spectrum in the MAGIC range to lower energies is also 
shown indicating that: i) there is a potentially smooth connection between the Fermi /LAT 
and MAGIC extrapolated data in the 3 to 10 GeV region, ii) the photon index steepens 
from 1.9 in the Fermi /LAT range to 2.7 in the MAGIC range. These results agree with the 
analysis of wider temporal intervals during this flare and during the whole active period, in 
which the source spectrum is well described by a broken power law with an energy break 



falling between 1 and 3 GeV (Tanaka et al. 



20111 1. Furthermore it is found that the high 



energy tail (E> 2 GeV) of the Fermi /LAT spectrum of PKS 1222+21 extends up to 50 
GeV, with a photon index in the range 2.4-2.8. 
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vThz) 




logE(GeV) 



Fig. 3.— High energy SED of PKS 1222+21 during the flare of 2010 June 17 (MJD 55364.9), 
showing Fermi /LAT (squares) and MAGIC (circles) differential fluxes. A red bow tie in the 
MeV/GeV range represents the uncertainty of the likelihood fit to the Fermi /LAT data. 
The unfolded and deabsorbed spectral fit of the MAGIC data is also s hown as a r e d bow 
tie, extrapolated to lower and higher energies (dotted lines) according to lAbdo et al.l (120091). 
A thick solid line (photon index F = 2.7) indicates a possible extrapolation of the MAGIC 
deabsorbed data to lower energies. The thick dashed line represents the E BL absorbed spec 



trum o btained from the extrapolated intrinsic spectrum using the model by 
J2OI1L 



Dominguez et al. 



- 15 - 



4. Discussion 



4.1. EBL limits 



The interaction of very high energy 7-rays with low energy photons of the isotropic 
EBL is a process with an energy dependent threshold, thus leading to an imprint of the EB L 



density on the measured VHE 7-ray spectra of extragalactic sources ( iMazin &: Raud 120071 ) . 
For PKS 1222+21 {z = 0.432), the measured spectrum spans from 70 GeV to 400 GeV 
probing EBL photons in the range 0.1-1 /xm (i.e. UV to near infrared). 

The EBL constraints us ing VHE 7-rays are us ually derived assuming an intrinsic 



spectrum of the source (e.g. 



Aharonian et al. 



20061 ). In FSRQs, the presence of dense 
radiation fields of soft photons ca n lead to the internal absor ption of VHE 7-rays, mimicking 



harder-than-intrinsic spectra (e.g. 



Sitarek fc Bednarek 



20081 ). However, for realistic spectral 



distributions of the interri al photon fields it should not change the EBL limits significantly 



(iTavecchio fc Mazin 



200% . 



In our case the simultaneous data from Fermi /LAT, which is free from internal or 
external absorptions, has b e en used to constrain the intrinsic photon index in VHE (e.g. 



GeorganoDoulos et al 



one utilized by 



2010 



Finke fc Razzaaue 



Georganopoulos et al. 



2OO9I ) . We adopt a method similar to the 



(I2OIOI ): the intrinsic spectrum in the VHE regime is 
assumed to follow the extrapolation of the Fermi /LAT above 3 GeV with a F = 2.4. This is 
a conservative assumption since in reality the spectrum could soften with increasing energy. 



The upper limit (95% C.L.) on the optical depth, r,„ax, for VHE 7-rays can be obtained 



from: 



T,na>c(-E) = log 



_F^UE) - 1.64 • AF{E) 
where Fi,^t,{E) is the maximum intrinsic flux at energy E, F^^^^lE) and AF{E) are the 
MAGIC measured flux and its error, respectively. The maximum intrinsic flux has been 



(2) 
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normalized at 70 Ge V assuming the EBL model giving a maximum flux absorption of 30% 



flAlbert et al. 



2008a] ). The derived limits on the optical depth are shown in Fig. H] together 



with a compilation of the predicted optical depths for a source at z = 0.432 computed 
according to recent EBL models. The limits confirm previous c onstraints on the EBL 



models in the UV to near infrared regim es derived using VHE (|Aharonian et al. 



Mazin fc Raue 



2007 



Albert et al. 



T 



2006 



2008al ) and HE spectra (lAbdo et al.li2010aj ). Given the 
fact that the EBL models predict for this redshift a stronger absorption with increasing 
energy, our data do not indicate a softening of the spectrum within the energy range of our 



observations. 



4.2. VHE 7-ray emission 



In the framework of the currently accepted EBL models, the observed simultaneous 
VHE and GeV spectra are consistent with a single power law with index ~ 2.7 ± 0.3 
between 3 GeV to 400 GeV, without a strong intrinsic cutoff. This evidence suggests that 
the 100 MeV - 400 GeV emission belongs to a unique component, peaking at 2 — 3 GeV, 
produced in a single region of the jet. If the emission process is inverse Compton scattering 
on external photons by relativistic electrons in the jet, as commonly assumed, a strong 
softening of the spectrum is expected above few tens of GeV if the external photons derive 
from the BLR. This is due to the combination of two effects: the decreased efficiency of the 



IC scattering occuring in the Klein-Nishina (KN) regi me (e.g. 



and the absorption of 7-rays through pair production (IReimei 



2009 



Liu fc Bai 



20061) 



Ghise 



2007 



1 ni fc Tavecchio 



20m 



Tavecchio fc Mazin 



^We note that this absorption has been invoked by iPoutanen fc Stern! (120 lOf ) to explain 



the existence of an apparently "universal" break energy in the 7-ray spectrum of FSRQs at 
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Fig. 4.— Optical depth along the line of sight to PKS 1222+21 (z=0.432) for a range of EBL 
models and the limits (95% CL.) from the MAGIC measurement, assuming the limiting 



intrinsic photon index Fvhe 
determination as derived in 



2.4. The grey-s. 



laded area shows the uncertainties in the EBL 



Dominguez et al. 



(120111 ). sec. 6.1 and Fig. 13. 
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The energy above which the KN effects become important can be roughly expressed as: 
Ekn ^ 22.5 GeV, where 1^0,15 is the frequency of the target photons in units of 10^^ Hz 
(or E'kn — 75 A^m GeV in wavelength fim units). 7-ray absorption becomes effective when 
E^^ ~ 60u-l^ GeV {E^^ 200 A^^ GeV). Above this energy a cut-off is then expected. 
The importance of both effects in the 10-100 GeV band is reduced if the external photon 



field is associated with the 



scenarios (e.g. 



Sikora et al. 



R to rus (z/q = 10^^ Hz), as envisioned by the "far dissipation" 



20081 ). In that case both effects start to be important above 



1 TeV. The absence of a spectral break or cutoff in the spectrum observed by MAGIC, 
strongly suggests that the 7— ray emission is not produced within the BLR. 

The other important result of the MAGIC observation is the evidence of fast variability, 
tvar ~ 10 minutes, indicating an extremely compact emission region with transverse 
dimensions, i? ~ 1.3 x 10^''(5/10)(tvar/10 minutes) cm. This seems to be difficult to reconcile 
with the "far dissip ation" scenarios if the emission takes place in the entire cross section 



of the jet (see also 



Tavecchio et al. 



20101) ■ Estimatin g the size of the BLR, -Rblr from 
the accretion disk luminosity, L^isk = 5 x lO^'^ergs"^ (iFan et al.ll2006l ). the distance of the 
emitting region is expected to be around d > Rblr = 3 x 10^^ cm. Assuming a conical 



Marscher 



jet w ith constant opening angle 6^ (see however, the suggestion of recoUimation, 
19801 ). its size would be i? ~ 6jd ~ 3 x lO^^(0j/5°) cm. The absence of absorption features 
in the VHE spectrum allows also to exclude absorption within the emitting region and, 
together with the observed variability, t o put a lower limit to the Doppler factor of the 



source. From 



Dondi fc Ghisellini 



(119951 ). Equation (3.7), assuming a power-law photon 
index 1.5 for the spectrum of the optical target photons, we get a lower limit 6 > 15, in 
agreement with Doppler factors derived from radio observations (Section 1). 



A possibility to reconcile the spectral information (pointing to emission beyond 



2 GeV. 
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the BLR) and the fast variabihty is to invoke the presence of very compact emission 
regions embedded within the large scale jet, as already proposed by several authors to 



explain the exceptiona. 



( iGhisellini fc Tavecchio 



l y rap i d variability in 



2008 



Giannios et al. 



KS 2 155-304, MknSOl and AO 0235+164 



2009 



Marscher fc Jorstad 



20101). An 



alternative possibilit y is that the iet experienc es a strong recollimation forming a small 



304 (e.g. 



nozzle (e.g. 


Nalewaiko & Sikora 


2009 


) as already 


Brombere & Levinson 


2OO9I; 


Stawarz et al. 


2006) 



20061 ). Alternative scenarios involving 



proton-driven cascades or proton-synch rotron emission in 



generated by filamentation instabilities (iFrederiksen et al. 



amp 



ified magnetic fields, e.g. 



20101), could also play a role. 



In conclusion the MAGIC observations of VHE emission from the FSRQ PKS 1222+21 
put severe constraints on emission models of blazar jets. These results were obtained from 
a short observation of a flaring source thanks to the collaboration between the MAGIC and 
Fermi projects. Repeated and hopefully longer observations of flaring blazars with MAGIC 
and Fermi promise substantial progress in the study of extreme blazars. 
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